%% finite translation matrix
a=sym(diag(ones(1,9),1)+diag([1],-9));
for n=1:10
    disp('next n is');
    disp(n);
    a^n
end